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ABSTRACT 

The VERITAS collaboration reports the detection of very-high-energy 
(VHE) gamma-ray emission from the high-frequency-peaked BL Lac object 1ES 
1218+304 located at a redshift of z = 0.182. A gamma-ray signal was detected 
with a statistical significance of 10.4 standard deviations (10.4 a) for the ob- 
servations taken during the first three months of 2007, confirming the discovery 
of this object made by the MAGIC collaboration. The photon spectrum be- 
tween ~ 160 GeV and ~ 1.8 TeV is well described by a power law with an 
index of T = 3.08 ± 0.34 sta t ± 0.2 sys . The integral flux is $(E > 200 GeV) = 
(12.2 ± 2.6) x 10 _12 cm~ 2 s _1 , which corresponds to ~ 6% of that of the Crab 
Nebula. The light curve does not show any evidence for VHE flux variability. 
Using lower limits on the density of the extragalactic background light in the 
near to mid-infrared we are able to limit the range of intrinsic energy spectra for 
1ES 1218+304. We show that the intrinsic photon spectrum has an index that 
is harder than T = 2.32 ± 0.37 sta t- When including constraints from the spectra 
of 1ES 1101-232 and 1ES 0229+200, the spectrum of 1ES 1218+304 is likely to 
be harder than T = 1.86 ± 0.37 sta t- 

Subject headings: galaxies: active — galaxies: BL Lacertae objects: individual: 
1ES 1218+304 — gamma rays: observations 



1. Introduction 



One of the major discoveries of EGRET on the Compton Gamma Ray Observatory was 
the detectio n of high-energy emis sion from more than 60 active galactic nuclei (AGN) of the 
blazar class (IHartman et al.lll999l ). Blazars, which include BL Lac objects and flat-spectrum 
radio quasars (FSRQs), are characterized by non-thermal emission and their spectral energy 
distributions (SED) contain two broad peaks. The low-energy peak (radio to UV or X-rays) 
is commonly interpreted as synchrotron radiation from ultra-relativistic electrons moving 
along a plasma jet pointing towards the observer. The origin of the second peak is less 
certain. Several models, from pure leptonic or hadronic models to leptonic/ hadronic hybrid 
models, can explain the high-energy peak (X-rays to TeV gamma-rays) (see iBottcherl 120071 . 
and references therein). 

In the TeV energy range, 19 blazars and one radio galaxy (M 87) have be en established as 
emitt ers of TeV gamma-rays. High-frequency-peaked BL Lac (HBL) objects (jPadovani &: Giommi 
19951 ) are a subclass of blazars characterized by a synchrotron peak at X-ray energies, unlike 
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quasars that have higher luminosity and a synchrotron peak at optical/infrared energies. 
With the exc eption of the recen tly discovered low-frequency-peaked BL Lac (LBL) object 
BL Lacert ae (Albert et al.l 120071 ). the intermediate-frequency-peaked BL Lac (IBL) objec t 



W Comae (jSwordy 



20081 ). and the flat spectrum radio quasar 3C 279 ([Albert et al. 



2008b|) 



all TeV blazars are high-frequency-peaked BL Lac objects. Some of the se objects can show 



rapid, down to a few minutes, flux variabil ity at TeV energies (see e.g., iGaidos et al.lll996 
Aharonian el al.ll2007cl : [Albert et al1l2008ah . 



Broadband observations from the radio to very-high-energy (VHE) gamma-rays are nec- 
essary to understand the physics of the jets and emission mechanisms. VHE observations 
can also help constrain the intensity and spectrum of the extragalactic background light 
(EBL) (IStecker et al.lll992l ). which are important parameters for cosmologists to test our 
understanding of structure and star formation in the Universe. Produced by stars and par- 
tially reprocessed by dust, the EBL at near to mid- infrared ( near-IR, mid-IR) wavele ngths 
is a strong absorber of TeV gamma-rays via pair production (IGould fc Schrederlll967l ). Di- 
rect measurements of the EBL in this range are particularly difficult due to the dominant 
foreground of zodiacal light. The measurement of high-quality energy spectra for blazars 
over t he energy range from 100 GeV to 10 TeV can be used to gain information about the 
EBL ( IDwek fc Krennrichl 120051 ) . The analysis of TeV gamma-ray energy spectra of several 
blazars suggests that the intensity of the EBL in the near-IR to mid-IR band is close to the 
lower limit from source coun t s measured by the H ubble and Spitzer space telescopes (see 
Aharonian et alJbood . l2007bl : iMazin fe RauelboOTT ). 



Based on SED modeling and BeppoSAX X-ray spectra, several of the HBLs were pre- 
dicted to be TeV sources and s everal of them have indeed been detected at TeV energies 
JCostamante fc Ghisellin] 12002^ . 1ES 1218+304 is an X-ray-bright (F lkeV > 2/iJy) HBL 



object located at a redshift z = 0.182 (IBade et al.lll998l ). and it was predicted to be a good 



TeV candidate based on the position of its synchrotron peak at high energy and sufficient 



(Albert et al. 


2006). This object was 


in May 2006 ( 


Aharonian et al. 


2008a) 



responding to 
was reported. 



6 times the flux determined from an extrapolation of the MAGIC spectrum 



In this paper we report on the detection of 1ES 1218+304 in VHE gamma rays with 
VERITAS. A description of the telescopes and a summary of VERITAS observations of 1ES 
1218+304 are presented in §2. The results of the data analysis are presented in §3, and a 
discussion of our results in the context of the EBL density is presented in §4. 
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Observations 



The VERITAS observatory consists of an array of four 12-meter diameter imaging atmo- 
spheric Cherenkov telescopes (IACTs) located at the Fred Lawr ence Whipple Obser vatory 



(31°40'30" N, 110 o 57'07" W, 1268 m a.s. l) in southern Arizona flWeekes et all 120021 ) . The 



telescopes use the Davies-Cotton design (IDavies Cotton! 1 195 71 ) and each utilize 350 front- 
aluminized and anodized hexagonal glass facets with a total mirror area of 106 m 2 . Each 
camera consists of 499 photomultiplier tubes (PMTs) separated by 0.15° and covers a 3.5° 
field of view. Light concentrators reduce the dead space between PMTs and decrease the 
amount of ambient light seen by the PMTs. The analog signals from the PMTs are pre- 
amplified in the camera before being sent through coaxial cables to an electronics room 
located at the base of each telescope. VERITAS uses a three-level trigger system to reduce 
the rate of background events caused by fluctuations in the night sky light and by cosmic- 
ray showers while retaining multi-telescope images consistent with gamma-ray showers. The 
analog PMT signals are digitized using custom-designed 500 mega-sample-per-second flash- 
analog-to-digital converters (FADCs) and the data are archived to disk. Additional technical 



Holder et al. 


(2006 


) and 


Maier et al. 


(2007) 



The first two telescopes were operated in the stereoscopic observation mode from March 
2006 and the third and fourth telescopes came online in December 2006 and April 2007, 
respectively. VERITAS observed 1ES 1218+304 from January to March 2007 using three 
telescopes. The data were taken in wobble mode where the source is offset from the center of 
the field of view by 0.5° and the background is measured directly from different regions in the 
same field of view but away from the source region. After removing data taken under poor 
sky conditions or affected by various detector problems, we were left with a total observation 
time of 17.4 hours covering a range in zenith angle from 2° to 35°, with an average zenith 
angle of 14°. 



3. Data Analysis & Results 



The analysis of the data was performed using independent analysis packages (see lDaniel et al 

20071 . for details on the analyses). All o f these analy ses yield consistent results. After cal- 
culating the standard Hillas parameters (jHillaslll985l ). images with an integrated charge less 
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than 400 digital counts^ or with a distance from the center of the camera larger than 1.2' 
were rejected. The location of the shower direction in the field of view a nd the impact param 



eter of the shower core w ere calculated using stereoscopic techniques (jfiofmann et al.lll999 



Krawczynski et al.l 120061 ) . Events originating from a circular region with radius 9 = 0.158° 
centered on the position of 1ES 1218+304 are taken as the source region, and the background 
was estimated from the same field of view using the reflected-region and ring background 
methods, as described in iBerge et al.l (120071 ). A set of scaled cuts on the width and length 
parameters were used to identify gamma-ray events in the data (jKonopelko et al.lll995l ). The 
cuts were optimized a priori with Monte Carlo simulations of gamma-ray and hadron induced 
air showers. The cuts used here are (-1.25 < mean-scaled width/length < 0.5). An extensive 
set of Monte Carlo simulations was used to generate lookup tables to calculate the energy 
of the primary gamma rays. The effective area of the detector as a function of zenith angle 
and gamma-ray energy was also calculated from these simulations. The variability of the at- 
mospheric conditions and the overall photon collection efficiency are the major contributors 
to the systematic error of the energy estimation. 

Figure [T] shows the number of excess events calculated using the reflected-region method 
as a function of the squared angular distance (6* 2 ) between the reconstructed shower directions 
and the nominal position of 1ES 1218+304. A clear excess is visible below the angular cut of 
0.158°, correspond i ng to a statistical significance of 10.4 standard deviations using equation 
(17) in iLi and Mai (119831 ) (617 signal events, 1466 background events with a normalization 
of 0.25). Figure [2] shows a two-dimensional sky map of significances for a region centered o n 
the radio coordinates of 1ES 1218+304 (12 h 21 m 21.9 s , +30°10'37" J2000) faecker et al.lll995h . 
The excess is compatible with that expected from a point source. Fitting a two-dimensional 
normal distribution to the uncorrelated excess map yields a peak position with coordinates 
(12 /l 21 m 26.3 s ±2.8f iat ±7.0fj /s , +30°li , 29"±33 / ; iat ±l , 30'4 s J2000), in good agreement with the 
radio coordinates. The systematic uncertainty on the tracking of the VERITAS telescopes 
was ±90 arcseconds at the time the data were taken on 1ES 1218+30fl 

Figure [3] shows the light curve of the integral flux above 200 GeV for the months of 
January, February and March 2007. The average integral flux was calculated for each day 
assuming a spectral shape of dN/dE oc E~ r with F = 3.08. A statistical test (x 2 /dof = 
6.7/11) indicates that no statistically significant variability was detected in the data. 



1 An integrated charge of 400 digital counts calculated using a 10 ns integration window corresponds to 
~75 photoelectrons. 

2 Pointing monitors were installed on all 4 telescopes during the summer of 2007. This is expected to 
improve the pointing accuracy significantly in the future. 
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Figure H] shows the time-average differential energy spectrum for gamma-ray energies 
between 160 GeV and 1.8 TeV. The shape is consistent with a power law (x 2 /dof = 2.1/5) 
dN/dE = C x (E/0.5TeV)~ r with a photon index T = 3.08 ± 0M stat ± 0.2 sys and a flux 



normalization constant C = (7.5 ± l.l stat ± 1.5 



sys, 



x 10 



-12 cm 2 g 



L TeV 1 . The integral flux 



is $(E > 200 GeV) = (12.2 ± 2.6) x lO^cm-V 1 which corresponds to ~ 6% of the flux 
of the Crab Nebula above the same threshold. The MAGIC collaboration measured the 
differential energy spectrum of 1E S 1218+304 betwee n 87 GeV and 630 GeV during 6 nights 
of observations in 2005 January ([Albert et al.ll2006l ). No evidence for time variability was 
detected. The VERITAS measurement extends the differential energy spectrum to 1.8 TeV 
and is statistically consistent with the results reported by the MAGIC collaboration. 



4. Discussion and conclusions 



Due to the relatively large redshift of 1ES 1218+304 (z = 0.182) we expect significant 
attenuation from the interaction of high-energy gamma-rays with the low-energy photons of 
the EBL. TeV gamma-rays from extragalactic sources are expected to interact with optical 
low energy photons as they travel through inter-galactic space, leading to a cutoff in the 
measured gamma-ray spectrum. The optical depth r of the attenuation is a complicated 
function of the gamma-ray photon energy, the distance to the source (redshift z), and the 
cross-section for pair production, and it is related to the density and spectral energy distri- 
bution of the cosmic background radiation. The physical importance of the EBL lies in its 
relation to galaxy formation and evolution and to the star formation history of the Universe. 
If the intrinsic spectra of blazars extend to ~ 10 TeV, absorption features in the blazar 
spectra can be used to learn about the EBL in the mid-IR band. TeV spectra of blazars 
therefore have the pote ntial of pr o viding independent constrai n ts on the infrared background 
light density (see, e.g., iPrimackl (120011 ); ICoppi fc Aharonianl ( 119991 )). In addition to inter- 
galactic absorption, the measured blazar spectra are shaped by intri nsic absorption in th e 
blazar, as well as by Compton scattering in the Klein-Nishina regime (IModerski et al.ll2005l ). 



Intrinsic spectra of blazars are not known a priori. Nevertheless, it is possible to use 
the EBL lower limits from galaxy counts to restrict the range of EBL scenarios/models 
and discern a corresponding range of intrinsic blazar spectra. We present possible intrinsic 
energy spectra of 1ES 1218+304 that are compatible with EBL lower limits from galaxy 
counts. Furthermore, we broaden our study and include the energy spectra of 1ES 1101-232 



and 1ES 0229+200 as measured by H.E.S.S. (see lAharonian et al.ll2006l . l2007bl . respectively) 
to reach a conclusion about the intrinsic energy spectra of this small sample of TeV blazars. 
The combination of energy spectra in the sub- TeV to TeV waveband at large redshifts 
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(1ES 1101-232 and 1ES 1218+304 are two of the most distant HBLs) with energy spectra in 
the multi-TeV regime at a moderate redshift (1ES 0229+200) provides additional sensitivity 
to EBL spectra and the relative intensities in the near-IR and mid-IR. Consequently, this 
also leads to stronger constraints for the intrinsic spectra of these blazars. The combination 
of these three blazars allows us to make a statement about the hardness of blazar spectra in 
general, and sheds light on the question: how hard are intrinsic TeV blazar spectra? 

In order to obtain the intrinsic spectrum of a blazar, the measured spectrum must be 
corrected by unfolding the effects of the EBL. Figure shows the range of EBL scenarios con- 
sidered for unfolding the intrinsic spectrum of 1 ES 1218+304. The lowe r limits from galaxy 
counts derived from the HST deep sky survey (IMadau &: Pozetti 1 120001 ) provide the lowest 
possible EBL in the optical to near-IR (filled triangles). Galaxy counts from the Spitzer 
infrared observatory (IFazio et al.l 12004 ) in the mid-IR are complement ary (open quadran- 
gles) . A convenient parameterization of EBL scenarios, provided by iDwek fc Krennrich 
(120051 ). is used in this study for providing a limit to the hardness of the blazar spectrum of 
1ES 1218+304 as it provides a wide range of EBL spectra with different near-IR to mid-IR 
ratios. The EBL scenarios are parameterized using polynomials and are constrained by de- 
tections and lower limits in the UV, optical, and sub-millimeter wavebands. Further details 



can be found in IDwek fc Krennrich! (120051 ). which presents 12 possible EBL realizations. 
The optical depth was calculated using a flat universe cosmology with H = 70 km s _1 , 
Q m = 0.3, and f^A = 0.7. For illustrative purp oses we show two low EBL scenarios in Figure 
[S] parameterized by lDwek &: Krennrich! (120051 ) . The first, called LLL (dotted line) is an EBL 
realization representing a low-intensity stellar component, a low-intensity 15/xm EBL flux, 
and a low-intensity of the far-IR flux, and falls significantly below the galaxy counts in the 
mid-IR. The second, called LHL0.70 (dashed line) is a scaled version of the LHL realization 
and falls below galaxy counts in the optical to near-IR. Furthermore, we have considered a 
large variety of EBL scenarios with different spectral indices and shapes that are within the 
boundaries of the shaded area in Figure The upper bound is somehat arbitrary and was 
chosen to be the HHH (High near-IR, mid-IR and far-IR) scenario. 

Furthermore, for reference, we also derive the absorption-corrected intrinsic gamma-ray 
spectra for theoretical EBL models by Stecker, Malkan, and Scully (2006, hereafter SMS). In 
these models, the intergalactic IR photon flux and density were calculated using a backward- 
evolution method which started with existing galaxy populations and modeled the luminosity 
evolution of these galaxies back in time. It is, however, important to note that in this paper 
we do not adopt any particular model, we simply use constraints from EBL measurements 
to derive a lower limit to the hardness of the 1ES 1218+304 spectrum. 



Table 1 shows the intrinsic energy spectra photon indices of 1ES 1218+304 for a range 
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of EBL scenarios. The ones that fall below the limit from galaxy counts (see Figure [5]) are 
marked with an asterisk in Table 1 and are not viable, e.g., the LHL0.70, the LLL and the 
LLH scenarios. Scenarios that are still compatible with the lower limits from galaxy counts 
represent the softest possible gamma-ray spectra. When considering 1ES 1218+304 by itself, 
the softest intrinsic spectrum is described by a power law with dN/dE oc E _2 ' 32±a37stat and 
is derived from a scaled version of the LHL scenario (LHL0.76). 

However, when applying this analysis to previously reported blazar spectra from 1ES 1101- 
232 and 1ES 0229+200, the LHL0.76 would require an extremely hard intrinsic spectrum 
for 1ES 0229+20C§. Therefore, in order to provide a limit to the hardness of the blazar 
spectra based on these three sources we search for the softest possible intrinsic spectrum. 
As can be seen from Table 1, a search for the softest possible blazar spectrum among this 
sample of three blazars yields an EBL scenario (AHA0.45) that still requires the spectrum of 
1ES 1101-232 to have a power law described by (dN/dE oc E~ L78±0 - 20stat ) and for the spec- 
trum of 1ES 1218+304 to be as hard as (dN/dE oc E- L86±a37stat ). All other EBL scenarios 
yield harder spectra for one of th ese three blazars. A detai led analysis of the VHE spectra of 
1ES 0229+200 was carried out by lAharonian et al.l (]2007bl ). This analysis supported an EBL 
spectrum and density close to the lower limits from the Spitzer measurements, and a hard 
intrinsic spectrum for the blazar. Figure [6] shows the measured spectrum of 1ES 1218+304 
and two possible intrinsic blazar spectra that are compatible with the limits from galaxy 
counts. 

These results clearly indicate that blazar spectra, as evidenced by two of the most 
distant blazars (1ES 1101-232 and 1ES 1218+304), are hard. These va lues are still w i thin t he 
acceptable range predicted for shock acceleration in blazars (see, e.g.. IStecker et al.l (120071 )). 
however they are also close to a limit (dN/dE oc E~ L5 ) that was previously suggested by 
Aharonian et al.l (120061 ) . The hard intrinsic spectrum of 1ES 1218+304 indicates that the 
peak in the VHE power output is located beyond ~ 2 TeV. A high-energy peak above 
2 TeV in the blazar spectral energy distribution was simi larly seen by H.E.S.S. f or both 
1ES 1101-232 (kharonian et alil2007ah and 1ES 0229+200 (| Aharonian et alil2007bh . 



In conclusion, results presented here fro m VERITAS observ ations confirm with high 
statistical significance the MAGIC discovery (lAlbert et al.l 120061 ) of the HBL object 1ES 
1218+304 as a source of VHE gamma rays. The normalization of the flux and spectral 
index both agree within errors with the MAGIC results. New limits on the density of the 



3 The data points for the energy spectra of 1ES 1101-232 and 1ES 0229+200 were provided to us by the 
H.E.S.S. collaboration. We have carefully checked the results of our power law fits against those published 
by the H.E.S.S. collaboration and they are in agreement. The correction for the EBL absorption was applied 
to the individual flux points and we fitted the absorption-corrected flux points with power laws. 
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EBL in the near-IR to mid-IR could not be established. However, based on lower EBL 
limits from galaxy counts we were able to limit the range of intrinsic energy spectra for 
1ES 1218+304 and showed that the intrinsic spectrum is harder than a power law with 
dN/dE oc E~ 2 - 32±a37stat . When including constraints from the spectra of 1ES 1101-232 and 
1ES 0229+200, the spectrum of 1ES 1218+304 is inferred to be harder than dN/dE oc 
g-i.86±o.37 Bta t _ Future deep observations at large zenith angles (where the effective area 
is larger for multi-TeV gamma-rays) could help extend the spectrum to 10 TeV and help 
constrain the EBL. We see no evidence of a high energy peak in the SED for 1ES 1218+304 
up to ~ 2 TeV; simultaneous, broadband measurements of the SED are required for detailed 
blazar modelling studies and to distinguish between leptonic and hadronic blazar models. 
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Table 1. The absorption-corrected spectral indices (dN/dE oc E~ rsourcc ) of 1ES 1101-232 

and 1ES 1218+304 and 1ES 0229+200 from a small sample of all EBL scenarios 
considered. EBL scenarios that fall below the lower limits from galaxy c ounts are marked 



with an asteri sk. Sc aled Aharonian scenarios (AHAx.y) are taken from lAharonian et al. 



(120061 ). The errors represent only the statistical uncertainty. 



Scenario 




r 

L source 




1ES 1101-232 


1ES 1218+304 


1ES 0229+200 


AHA1.0 


0.44 + 0.20 


0.35 + 0.40 


2.35 + 0.14 


AHA0.55 


1.54 + 0.20 


1.59 + 0.37 


2.41 + 0.14 


AHA0.45 


1.78 + 0.20 


1.86 + 0.37 


2.43 + 0.13 


HHH 


-0.67 + 0.12 


-0.72 + 0.29 


0.90 + 0.17 


LLH* 


2.01 + 0.22 


2.07 + 0.35 


2.12 + 0.20 


LHL 


2.04 + 0.20 


2.08 + 0.39 


0.94 + 0.32 


LHL0.70* 


2.32 + 0.21 


2.43 + 0.37 


1.43 + 0.29 


LHL0.76 


2.23 + 0.21 


2.32 + 0.37 


1.30 + 0.29 


LHL0.82 


2.18 + 0.21 


2.26 + 0.38 


1.20 + 0.30 


MHL0.70 


1.26 + 0.19 


1.34 + 0.36 


1.35 + 0.21 


MHL0.55 


1.61 + 0.19 


1.73 + 0.35 


1.59 + 0.20 


LLL* 


2.06 + 0.16 


2.20 + 0.34 


2.11 + 0.20 


SMS baseline 




1.70 + 0.34 
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Fig. 1. — Distribution of excess events as a function of the squared angular distance (@ 2 ) 
between the reconstructed shower directions and the nominal position of 1ES 1218+304. The 
dashed line indicates the boundary of the signal region at 0.158° (0.025 deg 2 ) as determined 
from optimization on the Crab Nebula. The solid curve shows the expectation for a point 
source as measured from data on the Crab Nebula. 
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Fig. 2. — Significance map of the region around 1ES 1218+304. The white cross indicates the 
position of the corresponding radio source. The white circle in the upper left corner shows 
the angular resolution of VERITAS (PSF). The ring background model w as used to e s timat e 
the background and the significances were calculated using equation 17 in lLi and Mai (119831 ). 



-16- 



CM 



V> 40 



E 

o 



CM 



S 30 



> 
o 

o 
o 

CM 
A 
LU 

o 



20 



10 



■10 



J I L 



J I L 



J I L 



54120 



54140 



54160 54180 
time [MJD] 



Fig. 3. — Light curve of the integral photon flux above 200 GeV for the source 1ES 1218+304. 
Each point corresponds to the average daily flux and assumes a spectral shape of dN/dE oc 
E~ r with T = 3.08. The error bars represent the statistical uncertainty. 
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Fig. 4. — Differential energy spectrum of VHE photons above 160 GeV for 1ES 1218+304. 
The markers indicate measured data points and the continuous line is a p ower-law fit. D own- 
ward pointing arrows correspond to upper flux limits (99% probability, Helene ( 19831 )) for 
bins with significances below two standard deviations. The dot - dash line shows t he differ- 
ential energy spectrum of 1ES 1218+304 measured by MAGIC (lAlbert et al.l 120061 ) and the 
dash line shows the differential energy spectrum of the Crab Nebula measured by VERITAS 
(September to November 2006) for comparison. 
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Fig. 5. — The range of EBL densities that are considered for reconstructing the source 
spectrum o f 1ES 1218+304 is shown . The filled triangles represent lower limits from galaxy 
counts by (IMadau fc Pozetti 1 120001) w hereas the open rectangles are lower limits from the 



Spitzer observatory (IFazio et al. 



200J). 
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Fig. 6. — Measured spectrum of 1ES1218+304 (circles) and de-absorbed spectrum for the 
AH AO. 45 EBL scenario that produces the softest possible spectrum among all EBL scenarios 
considered (filled quadrangles). Furthermo r e we show the correc ted spectrum for the SMS 
baseline model by lMalkan fc Steckerl (I1998I ); IStecker et all (120061 ) (hollow quadrangles). 



